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ABSTRACT 

One of the challenges to increasing the mass of a white dwarf through accretion is the tendency for the 
accumulating hydrogen to ignite unstably and potentially trigger mass loss. It has been known for many years 
that there is a narrow range of accretion rates for which the hydrogen can burn stably, allowing for the white 
dwarf mass to increase as a pure helium layer accumulates. We first review the physics of stable burning, 
providing a clear explanation for why radiation pressure stabilization leads to a narrow range of accretion rates 
for stable burning near the Eddington limit, confirming the recent work of Nomoto and collaborators. We 
also explore the possibility of stabilization due to a high luminosity from beneath the burning layer We then 
examine the impact of the /5-decay-limited "hot" CNO cycle on the stability of burning. Though this plays 
a significant role for accreting neutron stars, we find that for accreting white dwarfs, it can only increase the 
range of stably-burning accretion rates for metallicities < O.OIZ©. 

Subject headings: accretion, accretion disks — binaries: close — instabilities — novae, cataclysmic variables 
— nuclear reactions, nucleosynthesis, abundances — stars: white dwarfs 



1. INTRODUCTION 

The commonly accepted explanation for Type la super- 
novae is the thermonuclear ignition of carbon in the core of a 
massive white dwarf (see review of iHillebrandt & Niemeveii 
I2OOOI) . The favored channel for mass growth of the 
white dwarf (WD) is accretion from a companion star, a 
topic that has been well-studied (e.g., Sienkiewicz 1975; 
iPaczvnski & Zvtkowlll978l: ISion et al.lll979t INomotolll982l) . 
However, thermal instabilities in the accreting hydro- 
gen/helium layer result in thermonuclear runaways that 
may inhibit mass growth for accretion rates, M, below 
a critical value of ~ 1 0"^ M( 7) yr~' (jP aczvnski & Zytko' 
i978; Sion et alj ll979l: ISi enkiewk3 IT980I: lFuiimotoHl982 
Paczvnski 19831 : lYoon et al...2004,) . It is clear that for much 



bility criterion for hydrogen accretion onto a WD and pro- 
viding an explanation as to why the calculated range of M's 
for stable burning is so narrow. In ^ we incorporate the full 
CNO cycle, including the temperature-insensitive /9-decays, 
and show that this can modify the boundary for stable burn- 
ing only at low metallicities Z < O.OIZ©. These metallicities 
are so low that it is unlikely they are relevant to most progen- 
itor scenarios. We close in |5]by summarizing our work. 

2. STEADY-STATE HYDROGEN-BURNING 
SOLUTIONS 
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lower rates < 10 ' Mq yr ^ the accumulation of the hydro- 
gen layer always leads to mass ejection during a classical nova 
event (F uiimoto 1982a b; Macdonald 1983; Gehrz et al. 1998|; 
iTownslev & Bildstenll2004 ). 

Analytic studies of steadily burning and accreting atmo- 
spheres on WDs have neglected the /3-limited CNO cycle, 
whereas studies of neutron stars find t hat the CN O cycle is 
strong ly limited by /3-decays (Fuiimoto et al.|[T9 81; Bildsten 
1 19981) . When the rate-limiting step is a /?-decay, the nuclear 
energy generation rate is independent of temperature, and thus 
the burning is stable to thermal perturbations. What is not 
known is the metallicity of the material accreted on the WD 
needed for the /3-limited CNO cycle to become important dur- 
ing stable burning. If the inclusion of the full CNO cycle can 
meaningfully increase the range of accretion rates that lead to 
thermally stable burning layers on accreting WDs, it would 
substantially impact our understanding of Type la progenitor 
scenarios. 

We begin in ^and ^by rederiving the thermonuclear sta- 



We follow the one-zone formalism of iPaczvnskil ( Il983h . 
in which quantities such as temperature and pressure have a 
single value and are defined at the bottom of the hydrogen- 
burning layer. We neglect the effect of the unde rlying he- 
lium la yer except as a source of luminosity, but see lJose et alj 
(11993b for a numerical two-zone study of unstable hydro- 
gen and helium shell flashes on WDs, and lCooper & NaravanI 
(.2006) for a two-zone model of X-ray bursts on neutron stars. 
We assume a non-degenerate equation of state, and we include 
radiation pressure. We also assume radiative energy transport 
via Thomson scattering. These assumptions are justified in 
Figure [1] 

By energy conservation in steady-state, the exiting radiative 
luminosity, Lq (the 0-subscript denotes the steady-state value), 
must be equal to the sum of the luminosity emitted from the 
core, L^, plus the luminosity from nuclear burning. It is safe 
to ignore the luminosity from compression of the accreting 
material as it is or ders of magni tude lower than the luminosity 
from burning (see lTownslev & B ildsten (2004) for an analytic 
approximation to the compressional luminosity). In steady- 
state, material is burned as fast as it is accreted, so 



Lq = MXE+L^, 



(1) 
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where the hydrogen mass fraction is taken to be X = 0.7 unless 
otherwise specified, and E = 6.4 x 10'^ erg g"' is the energy 
released when a unit mass of hydrogen is converted to helium. 
The radiative-zero solution gives the luminosity as 



47r GMc aT^ 
1 P~' 



(2) 



where the Thomson opacity is kt = (1 +X)aT: /2m^ and M is 
the mass of the WD. The envelope mass, Menv, is very small 
and is neglected in comparison with M. 

The natural unit for the luminosity is the Eddington lumi- 
nosity. 



AttGMc 



(3) 



Since stable nuclear burning on a WD releases much more 
energy than gravitational infall, the accretion rate is given in 
units of the nuclear Eddington accretion rate assuming com- 
plete hydrogen depletion. 
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AttGMc 



(4) 



ktXE 

Equations ([1]) and (|2]l are rewritten in terms of /?, the ratio 
of gas pressure to total pressure, as 



l-/?o = fl + /, 



(5) 



where the dimensionless accretion rate and core luminosity 
are d = M /Meh and / = Lc/LEdd, respectively. Equation 
dH) is an extension of the Eddington standard model (e.g., 
[Hansen et al.ll2004 l). This relation is useful when we discuss 
the effect of radiation pressure on stability. Since < /3 < 1, 
a necessary condition for steady-state burning is d + f < 1 ; 
that is, the sum of the contributions to the luminosity from 
accretion and the core cannot exceed the Eddington limit. 
The specific entropy, s, is governed by 



ds 
T—-- 
dt 



dL_ 
dM 

47r/?2 P' 



(6) 



where e is the energy generation rate, and the burning layer 
mass has been approximated as the mass contained in a scale 
height, h = P/pg- The gravitational acceleration at the base 
of the envelope is g = GM/R^, with R the radius of the WD 
core. Using a standard thermodynamic relation, equation ^ 
is rewritten in the form 



dT 
dt 



1 

cp 



L-L^ g 
AttR^ P 



T dP 

+ -Vad— . 

P dt 



(7) 



The specific heat at constant pressure when non-degenerate is 

32-24/3-3/32 



^ 2/3^ /imp ' 

where /i is the mean molecular weight, and 



91nr 



'ad ■ 



ainP 



2(4-3/3) 
32-24/3-3/32 



(8) 



(9) 



is the adiabatic temperature gradient. The first term on the 
right side of equation (Q contains contributions from heat- 
ing by nuclear burning and core luminosity into the zone, and 
cooling by radiative luminosity out of the zone. The second 



term represents a temperature change due to adiabatic expan- 
sion or compression of the accreted matter. 

In steady-state, the time-derivatives in equation (|7| are 
equal to zero, so 



£0 = 



^Q-^c g 

AttR^ Po' 



(10) 



Thus, given Lc and M, we can find steady-state (but not nec- 
essarily stable) solutions for 7b and po, as shown in Figure [T] 
for the case of hydrogen-burning in the "cold" CNO limit. 

2.1. The Cold CNO Limit 

Hydrogen burning proceeds via the CNO cycle at the 
temperatures ( > 10^ K) relevant to accretion rates of M > 
10-'" Mq yr-> ( lFuiimotolll982bl:lTownslev & Bildstenll2004l) . 
The slowest step in the CNO cycle for T < 10^ K and p < lO"* 
g cm"^ is ''*N(p,7)'^0, and so nearly all catalytic nuclei are in 
the form of ''^N. In this "cold" CNO hmit, the energy genera- 
tion rate is 



Ccold = ■ 



6CN0 ^«14N 

p dt 



= 4.4 X lO^V^ZcNO X 



exp (-15.231/75 

7^2/3 



1/3 



I f 3.0057 A ' 

+8.3 X 10-^- " ^ ^ 1 — -' 



_exp| 



,3/2 



ergg s 



(11) 



where Qcno = 4.0 x 10"^ erg is the energy released for one 
loop of the CNO cycle, the reaction rate for '"^N is obtained 
from Lemut et al. (2006), and Tg = T/IO^ K. The mass frac- 
tion of CNO isotopes is Zcno, and we assume the total metal- 
licity to be Z = 3Zcno.' The first term in the parentheses is 
the non-resonant reaction with an S-factor of 1.7 keV-barn. 
The second represents the 259 keV resonant reaction, which 
dominates for T > 1.5 x 10^ K. 

Shown in Figure [1] are the locations of the steadily -burning 
hydrogen layer for WDs of mass 0.5 Mq (squares) and 
1 .35 Mq (circles), accreting matter at M > 10"^ M© yr"' with 
metallicities Z = lO'^ (solid lines) and Z = lO""* (dashed lines). 
The hydrogen is assumed to burn in the cold CNO limit. Filled 
points are thermally stable (see To the left of the long- 
dashed and long-dashed-short-dashed lines, the layer is non- 
degenerate and opacity is set by Thomson scattering, respec- 
tively. The steady-state solutions are well within this region, 
so our approximations are valid. Earlier ex amp les of these so- 
lutions can be found in lSienkiewiczl (119751) and lNomoto et alj 
( 12006.) . For completeness, we also show the location where 
steady-state helium burning ( ^3.51 ) would occur for a 0.9 Mq 
WD (pentagons). 

A few trends are clear For a fixed M, increasing the WD 
mass, and thus g, drives the steady-state temperature and den- 
sity higher The dotted line shows where the p+'^^N reaction 
timescale is equal to the '^O /3-decay time of 176 s; above 
this curve the CNO cycle is /^-limited. At low metallicities 
of \Q~^Zq (dashed lines), the layer is hot enough that the j3- 
limited CNO cycle is non-negligible for the highest accretion 

* Though there is evidence for enhanced CNO abundances in nova ejecta 
IGehrz et al.ll 19981) , we will presume that the metallicity in the accumulating 
layer matches that of the accreted material. 
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Fig. 1. — Temperature and density at the burning depth for varying WD 
mass and accretion rate assuming steady burning. Shown are WDs burning 
hydrogen in the cold CNO limit with M = 0.5 Mq (squares) and 1.35 Mq 
(circles), with metallicities Z= 10"^ (solid lines) and 10"'' (dashed lines); Z 
is not given in solar units. Also shown are solutions for a 0.9 Mq WD accret- 
ing pure helium (pentagons). The lower-right-most point of the hydrogen- 
burning solutions has M = 10"^ Mq yr"' , while the lower-right-most helium- 
burning solution has M= 10"^ Mq yr"'. Connected points moving towards 
the upper-left increase logarithmically with 8 points per decade of M. Filled 
points are thermally stable. Also shown are lines to the left of which the layer 
is non-degenerate (long-dashed line), opacity is due to Thomson scattering 
(long-dashed-short-dashed line), and radiation pressure dominates (dashed- 
dotted line), and below which hydrogen bums in the cold CNO limit (dotted 
line). 

rates and most massive WDs. Since these are the WDs of di- 
rect interest to the Type la progenitor problem, we pursue this 
question further in ^ However, we first wish to carry out 
the analysis needed to determine the thermal stability of the 
steadily-burning solutions. 

3. THERMAL STABILITY OF THE SHELL SOURCE 

We begin by perturbing the temperature around its steady- 
state value via 

T{t) = TQ + 6Tt'l'', (12) 

where 5T /To ^ 1 . If the perturbation timescale, |t|, is shorter 
than the accretion timescale, face = M^vi/M, then Menv can be 
approximated as constant during the perturbation. With this 
assumption, the accompanying pressure perturbation can be 
related to the density perturbation by the parametrization 

d\nP 



d\np 



(13) 



Mem 



as shown below. In our final results, the inequality |r| ^ face 
is indeed satisfied for essentially the entire parameter regime, 
so this assumption is self-consistent. 

3.1. Relations among Temperature, Pressure, and Density 
Perturbations for a Constant Mass Envelope 

For a constant luminosity envelope with Thomson opacity, 
P(r) (X T{r)^, so p{r) oc /"(r)^/"*, neglecting degeneracy. Since 
the envelope mass is much smaller than that of the WD, inte- 
gration of hydrostatic equilibrium yields the density profile 

3 



P(r) = Pb 



1-1 

X 



1-^ 
r 



(14) 



where x = 4h/ R\„ and the scale height is /i = Pb/ Pbg- (To avoid 
confusion, in this section we explicitly label variables at the 
bottom of the layer with the subscript "b".) The outer radius 
of the envelope, where p(Ro) = 0, is 



Ro = Rh/(l-x). 



(15) 



Thus, the scale height is constrained to be h < Rh/4-. The thin 
solid lines in Figures |2] and |4] are accretion rates above which 
this constraint is not met, i.e., where hydrostatic envelope so- 
lutions do not exist. This is a stronger bound than the nuclear 
Eddington limit (dotted line in Figure |2]l. This hydrostatic 
constraint is similar to the "red-giant lumin osity" dPaczvnskil 
Il970t lFujimo"tolll982bt iHachisu et alJI 19991) above which the 
envelope expands to red-giant dimensions and outfl ow may 
occur. The numerical results of 'Nomoto et alj ( l2006h for this 
upper bound are shown in Figure |2] (f/i'n dashed line). The 
discrepancy at lower WD masses is likely due in part to the 
approximation of the burning layer as the mass within a scale 
height. 
The mass of the envelope. 



47rr p(r)dr 



Rb 



3 1 

--4Tr pi,Rl^ 



In 



1- 



(16) 



is conserved during the perturbation, so (5Menv = 0. Thus, 



Slnx 



1- 



In 



1 



1- 



2 3 
X X 



= D. (17) 

From the definition of x, we have Slnx = SlnP^- Slnp^, so 
SlnPb = (1 + 1 /D)(51npb, and a = (1 + 1 /D). Thus, a is com- 
pletely specified by the steady-state burning conditions. Note 
that X < I constrains a to be < 1 as well. 

The equation of state yields the density perturbation in 
terms of the temperature perturbation: 

{dlnp/dlnT)p 



6\nph = - 



-S\nn. 



(18) 



l-a{dlnp/dlnP)^ 

For a non-degenerate gas, (91np/91nP)^ = 1//? and 
{dlnp/dlnT)p = 3-4/f3, so 



S\npi,=- 



SlnTb. 



(19) 



\l-a/(3 

3.2. Thermal Stability Condition 

We apply these perturbations to equation (|7|. The steady- 
state solutions will be stable if M[r] < 0. Keeping only first- 
order terms, we arrive at the condition for stability, 

A 



where 



B-- 



1 + Vaa(^ 



_ 4-3/3 
~ l-a/l3 
4-3/3 



(20) 



\-a/(3' 

The logarithmic derivatives of e are 
91ne 



A = 



d\np 



A a 



dint 



f 

2 + - 



dlnT 



(21) 



(22) 
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Fig. 2. — Range of accretion rates for which the hydrogen-burning layer 
is thermally stable, assuming Z = 10 - and no core luminosity. Above the 
thin solid line, hydrostatic envelope solutions do not exist. Below the thick 
solid line, the atmosphere is thermally unstable. Stable burning is only pos- 
sible between the two solid lines. Also shown are equivalent bounds (dashed 
lines) from a recent numerical simulation tNomoto et al..2Q06.) , as well as the 
nuclear Eddington limit {dotted line). 



Stability can be achieved via two physically distinct mech- 
anisms, when either A or B is negative. The case of A < 
(the thick solid line in Fig. |2]is A = 0) is the condition where 
an increase in temperature causes a larger loss of energy due 
to exiting radiation than gain from energy generation. This 
allows for stability when the shell is thin (h <C R). The 
second stability mechanism involves the parameter B, which 
is the ratio of gravothermal specific heat, c», to cp. Since 
cp > 0, the inequality B < is equivalent to c, < 0. Phys- 
ically, this is a situation in which the addition of heat to the 
shell causes a loss of internal energy due to expansion work 
dKippenhahn & Weigerllll990l) . Note that the layer is unstable 
if both A and B are negative, but, with the assumption that A 
and ly are non-negative, this case is ruled out. 

Increasing the fraction of radiation pressure (thereby de- 
creasing P) makes th e lay er more stable: this is radiation 
pressure stabiUzation dSugimoto & FuiimotollT978t iFuiimotol 
Il982bt I^onetan i2004r In broad terms, if the plasma is ra- 
diation pressure dominated, a small increase in temperature 
with constant or decreasing pressure will cause the density to 
decrease dramatically. This leads to a large drop in energy 
generation as well as a large amount of expansion work, both 
of which quench a thermal instability. From equation (|5]l, de- 
creasing /3 can be achieved by increasing the accretion rate 
and/or core luminosity. This helps to explain the requirement 
of high accretion rates and/or very hot WDs for stability. 

3.3. The Cold CNO Limit 

The thick solid line in Figure |2] is the minimum accre- 
tion rate (i.e., A = in eqn. ( l20l l) for which the hydrogen- 
burning layer is thermally stable as a function of M, assuming 
Z = 10"^ and no core luminosity. The cold CNO limit is not 
assumed, although it is still applicable. The slight upturn in 
the thick solid line for high masses is due to the 259 keV res- 
onance in the ''*N proton capture, which begins to become 
non-negligible near T = 1.5 x 10** K. The nuclear Eddington 
limit is shown for reference (dotted line). Nearly independent 



of mass, the lowest stable accretion rate is ~ 1/3 the limiting 
accretion rate above which no hydrostatic envelope solutions 
are found (thin soli d line). Also plotted a re equivalent bounds 
(dashed lines) from lNomoto et all ( l2006l) . We now derive the 
near constancy of this narrow range of thermally stable accre- 
tion rates. 

The upper bound (thin solid line) is where x = 1, which 
gives 



GM/R 



(23) 



where the subscript "max" refers to the quantity at the max- 
imum accretion rate. At the lower bound of stability (thick 
solid line), h <^R/A, and we expand a in powers of h/R to 
obtain 

Expanding equation dSTT i in terms of h/R, assuming cold CNO 
burning and no core luminosity, and only keeping the first or- 
der term, yields the condition at the lower bound. 



v=l + 



4 \6k^Ty^ia/ ump 



2+- 



1 



3 (25) 



Pmin 5 PtninGM/R \ Pmin J \ At 

where the "min" subscript refers to the value at the minimum 
accretion rate for stability. We have also set A = 1, which 
is appropriate for the cold CNO limit. From Figure [T] we 
see that the temperature does not change significantly over 
the range of stable accretion, so we approximate Tmm = Tinax- 
Then, equation ( l23T l gives 

4 /?„.x 1 \ / 4 



1 + 



4 ^ 

Solving for the ratio of accretion rates yields 



2+- 



1 



1- 



5/?2 



1-4/ 13, 



(26) 



(27) 



-/3min V 4 Z/ + 4-6/3mi 

Finally, we expand this ratio around /Jmi,, = 3/4 and, given that 
the logarithmic temperature derivative of the cold CNO cycle 
is 

5.077 2 
3 



,1/3 



(28) 



forr < 1.5 X 10* 

M,. 



K, we assume = 10 to obtain 



= 2.9 



1+2.0 /3, 



(29) 



The coefficient in front varies from 3.1-2.6 for v = 9-\2, and 
since P^i^ never falls below 1/2, the ratio of accretion rates 
is tightly constrained to be ~ 3, regardless of mass. This is 
the reason for the nearly constant separation between the two 
solid lines in Figure|2] 

3.4. Importance of Including the Pressure Perturbation 

To simplify the preceding analysis, one might be tempted 
to assume a constant pressure perturbation (a = 0) because 
the layer is "thin" - i.e., h <^R- for the minimum thermally 
stable accretion rate. However, this will result in significant 
errors, as we now demonstrate. 

Expanding equation ( 1211 1 in terms of a and only keeping the 
first order term yields the revised stability condition 



i4 



f 



< 0, (30) 
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Fig. 3. — Critical accretion rates vs. core luminosity for a 1.35M0 WD 
with the given metallicities (solid lines). The atmosphere is thermally stable 
to the right of these curves. Also plotted is the luminosity from steady burning 
of the helium below the hydrogen layer {dotted line). For a given accretion 
rate, this Une represents the maximum possible core luminosity for hydrogen 
accretion. 

or, with the assumption of cold CNO burning and negligible 
core luminosity. 



(31) 



Thus, even if the scale height is only 5% of the WD ra- 
dius for /3 2± 1/2, allowing a pressure perturbation changes 
the stability condition from v <9 io v < 10.5. Given the 
strong temperature dependence of the cold CNO cycle, this 
is a significant change in the minimum thermally stable ac- 
cretion rate. For example, a 1.0 Mq WD can accrete stably 
down to 1.7 X 10"^ Mq yr~', but if we incorrectly assume a 
constant pressure perturbation, the minimum accretion rate is 
3.0 X 10"^ Mq yr"', nearly a factor of 2 higher 

3.5. Pure Helium Accretion 

In the case of a WD accreting pure helium, the fuel burns 
via the triple-a reaction, given in Hansen et al. (2004) as 

5 1 X io**o^y-' 

e3a = t;:t^ exp(-4.4027/r9) erg g"' s"' , (32) 



7^3 

^9 



SO A = 2, and 



V3o 



4.4027 



(33) 



See Figure [U for steady-state conditions of a Q.9Mq helium- 
accreting WD {pentagons). For this mass, the range of stable 
accretion rates is 1.7-4.0 x 10"'' yr"'. 

3.6. Flux Stabilization 

If the WD interior is hot and emits a core luminosity com- 
parable to that from the nuclear energy release of the accreted 
material, the burning can be stabilized dFujimoto & Truranl 
Il982t lPaczynskil[T983h . The regions to the right of the soUd 
curves in Figure [3] are stabilized by this mechanism. For ex- 
ample, a 1.35 Mq WD with Lc = 10"''' erg s"' will stably burn 
accreting material at a rate of 10""^ Mq yr"', orders of magni- 
tude lower than the cases presented earlier As the core lumi- 
nosity increases, a larger range of accretion rates are stable. 



Our results differ from those of iPaczvnskil (Il983h due to his 
assumption of a thin shell (P = Msmg/^TTP?) and lack of in- 
clusion of the full CNO cycle, among other reasons. 

While a large core luminosity seems an attractive possibil- 
ity for stabilizing the burning layer, it is difficult to identify an 
energy source that would provide enough luminosity to affect 
stability during the accretion of > O.IMq typically needed to 
ignite the carbon in the core for a Type la supernova. The 
most obvious option, the luminosity from the cooling WD it- 
self, is orders of magnitude too small. Even steady burning 
of the helium below the hydrogen layer is inadequate to sig- 
nificantly increase the stable accretion rate range. Complete 
and stable burning of the helium gives a luminosity, L = MEue 
(dotted line in Fig. where = 5.8 x 10'^ erg g"'. For a 
1.35Mo WD, this additional luminosity only lowers the min- 
imum stable accret ion rate by 10%. 

Recent work by IStarrfield et all (120041) finds stable hydro- 
gen burning on a hot (L ~ 10^^ erg s"') accreting WD of 
M = 1 .35 Mq for accretion rates down to 1 .6 x 10"^ Mq yr"' . 
Figure|3]shows that this high luminosity is not responsible for 
stabilizing the hydrogen burning. The physics that stabilizes 
their models at such low accretion rate s remains unclear and is 
at odds with the more recent efforts of lNomoto et al.l ( 120061) . 

4. THE FULL CNO CYCLE 

If the temperature and density at the burning depth are very 
high, all fusion reactions in the CNO cycle occur more rapidly 
than the /3-decays. The energy generation rate is then deter- 
mined by the /3-decays of '"'O and '^O, which are independent 
of temperature and density, giving ehot = 5.8 x IO'^Zcno erg 
g"' s"'. (Even for the lowest metallicity cases we consider 
where Z = lO"*", p-p burning is still negligible for the relevant 
accretion rates.) In this /3-limited CNO cycle, hydrogen burn- 
ing is stable to temperature perturbations. We now derive the 
conditions for which the temperature-independent /3-decays 
of the full CNO cycle can stabilize the burning layer 

The energy generation rate of the full reaction rate network 

is 



(34) 



where / refers to the isotope, j differentiates between fusion 
and /3-decay reactions for a given isotope, and Qij is the en- 
ergy released in the specified reaction. In particular, the CNO 
cycle has a fork at '^N; when cold and/or diffuse, '^N nuclei 
preferentially undergo a /3-decay to '^C, and when hot and/or 
dense, they capture a proton to become '^O. The nuclear reac- 
tion timescales, f,j, are defined to satisfy the differential equa- 
tions 

dni ■sj—^ nil 
dt ~ ^ tk 

k 



(35) 



where k refers to the isotopes that can produce isotope /. With 
this definition of nuclear reaction timescales, the proton cap- 
ture timescales are given as 



1 



1 



capture " 



np(crv} Xp{av)' 



(36) 



The reaction rates, (ctv), are only fu nctions of tem perature and 
are obtained from NACRE (Angul o et all 1 19991) and LUNA 
(iLemut et alJl2006l) . 

If the thermal perturbations occur on timescales much 
longer than reaction timescales, the reaction rate network is 
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always in equilibrium for a given temperature and pressure. 
The time for a given nucleus to go once around the CNO cy- 
cle is 



^cycle 



CNO^CNO 

I4m„e 



(37) 



where we have chosen the mean molecular weight of the cat- 
alytic nuclei to be lAirip. We find that this timescale is always 
much shorter than the perturbation timescale for all metallic- 
ities and thermally stable accretion rates. Thus, drii/dt = 0. 
Then, at a given temperature and pressure, we can solve for 
all abundances in terms of one of the abundances, such as 
n\2c- Defining the weighting factor. 



rii/hj 



(38) 



we have 



nnclhic 

TiTn 



Defining another weighting factor, /i,(r, p) = ni/nnc, yields 



^ZcNo , > EijQufu(T,p) 

= XZcNoG(r, p). 



(40) 



Here, A, is the atomic mass of the isotope, and G(T,p) is a 
function that depends only on the temperature and density. 
This implies that u is also solely a function of temperature 
and density. 

With the analytic form of the energy generation rate, finding 
the steady-state solutions is as before. Specifying the accre- 
tion rate, the core luminosity, the CNO metallicity, and the 
WD mass fixes the temperature and density at steady-state. 
These in turn fix the value of Thus, given L^, Zcno, and M, 
we can analytically find values of M for which the atmosphere 
will satisfy equation ( |20] | and thus be thermally stable. 

Metallicity-dependent ranges of accretion rates for which 
the layer is thermally stable, with burning via the full CNO 
cycle, are shown in Figure |4] For accretion rates above the 
thin solid lines, hydrostatic envelope solutions do not exist, 
whereas below the thick solid lines, thermal perturbations in 
the atmosphere grow exponentially. Lowering the metallicity 
of the accreting material and increasing the WD mass increase 
the temperature and density of the burning layer and lead to 
shorter timescales for proton captures. This lowers the mini- 
mum accretion rate for thermal stability because i' decreases 
as the temperature-independent /3-decays of '"'O and '^O be- 
come more important to energy generation. This effect is also 
evident in Figure [3] The limiting case of hot CNO stabiliza- 
tion occurs in the hydrogen-burning shells of neutron stars, 
where burning proceeds via the /3-limited CNO cycle for solar 
metal licities at very sub-Eddington accretion rates (Bildsten 
119981) due to the high gravity. 

5. CONCLUSIONS 

We have confirmed that stable burning of hydrogen on an 
accreting WD requires accretion rates within a narrow range 
for solar metallicities (see Fig. |2]i and shown that this is a ro- 
bust result set by radiation pressure stabilization. Increasing 
the range of thermally stable accretion rates can be temporar- 
ily accomplished by having a luminosity exiting the WD core 
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Fig. 4. — Ranges of thermally stable accretion rates assuming no core 
luminosity, as in Figure |2] but with the given metallicities. Burning is via 
the full CNO cycle. The upper bound is the rate above which hydrostatic 
solutions do not exist (thin solid lines). The location of the upper bound has 
a weak dependence on metallicity; upper bounds for metallicities between 
solar and Z = 10"'' lie between the thin solid lines shown. The lower bound 
is the transition between thermally stable and unstable burning layers (thick 
solid lines). 



of a magnitude comparable to that released from nuclear burn- 
ing of hydrogen at that accretion rate (see Fig. O. However, 
there is no energy source available to the deep core of the WD 
that would allow for such stabilization during accretion of the 
0.1-0.2 Mq of material needed to reach core ignition for a 
Type la supernovae. Thus, the challenge of a narrow accre- 
tion rate regime for stable hydrogen burning of solar metallic- 
ity material remains. 

In the hopes of alleviating this bottleneck, we considered 
the stabilizing effect of /3-limited burning of hydrogen via 
the "hot" CNO cycle. This effect reduces the accretion 
rate needed for stable burning, but requires metallicities of 
< 10"^ Zq before it opens up much parameter space (see Fig. 
m. Such metallicities are relevant to a few globular clusters in 
our galaxy and o thers, and for t his re ason, calculations such 
as those done bv lNomoto et alj (l2006h should be extended to 
lower metallicities as they might prove relevant to observa- 
tions of accreting sources in globular clusters. 

However, while very low met allicities can be fou nd in these 
systems, we know jTremonti et a l. 2004; iLee et al.iii200 6) that 
the metallicity for most of the stellar mass in the nearby uni- 
verse is larger than 10"^ Zq, making it appear unlikely that 
low-metallicity environments will play a significant role in the 
overall Type la supernova rate. Full resolution of the impact 
of the hot CNO cycle awaits the combination of multi-zone 
calculations with a population synthesis code capable of as- 
sessing the contribution of stable burning at lower accretion 
rates to the Type la rate. 
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Nomoto, and Tony Piro for discussions. This research was 
supported by the National Science Foundation (NSF) through 
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